Optical absorption and energy-loss spectra of aligned carbon nanotubes 
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Optical-absorption cross-sections and energy-loss spectra of aligned multishell carbon nanotubes 
are investigated, on the basis of photonic band-structure calculations. A local graphite- like dielectric 
tensor is assigned to every point of the tubules, and the effective transverse dielectric function of the 
composite is computed by solving Maxwell's equations in media with tensor-like dielectric functions. 
A Maxwell-Garnett-like approach appropriate to the case of infinitely long anisotropic tubules is also 
developed. Our full calculations indicate that the experimentally measured macroscopic dielectric 
function of carbon nanotube materials is the result of a strong electromagnetic coupling between the 
tubes. An analysis of the electric-field pattern associated with this coupling is presented, showing 
that in the close-packed regime the incident radiation excites a very localized tangential surface 
plasmon. 
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I. INTRODUCTION 

The description of the electronic response of carbon 
nanotubesEfB has been a challenge for theoretical and 
experimental investigations. Various theoretical studies 
of plasmon excitations in single-shell carbon nanotubes 
were reported,trO electron-energy-loss spectra frCjift in- 
dividual multishell nanotubes were investigatedErEI by 
changing the number of shells, and curvature was found 
to induce little effect on the covalent bonding of multi- 
shell nanotubes W)A the availability of aligned car- 



bon nanotube films," H i optical measurements were car- 
ried out with polarized light ,E3ll3 thereby evaluating 
the frequency-dependent effective dielectric function of 
the composite and showing that carbon nanotubes have 
an intrinsic and anisotropic metallic behaviour. Accu- 
rate calculations of the effective dielectric function of 
densely packed carbon nanotubes, as obtained by solving 
Maxwell's equations with the use of tensor-like diele 
functions, have been carried out only very recently £i 

In this paper we report extensive calculations of the 
frequency-dependent effective dielectric function of a 
composite made up of aligned carbon nanotubes embed- 
ded in an otherwise homogeneous medium. In section 
II our effective medium theory is described. We take an 
electromagnetic wave normally incident on the structure, 
and focus on the case of electromagnetic waves polar- 
ized normal to the cylinders (p polarization). For this 
polarization w*i also introduce a generalized Maxwell- 
Garnett (MG)lia effective dielectric function appropriate 
for anisotropic tubules. Absorption of p polarized light 
is found to be sensitive to both the inner cavity of hol- 
low tubules and the anisotropy, which we first investigate 
in section IIIA within the MG approach. In section IIB 



we focus on the close-packed regime, showing that the 
experimentally measured macroscopic dielectric function 
of aligned carbon nanotubes is the result of a strong elec- 
tromagnetic coupling between the tubes. Calculations of 
the electric field, the induced charge, and the so-called 
energy-loss function, i.e., the imaginary part of the ef- 
fective inverse dielectric function, are also presented, for 
various values of the filling fraction and the ratio of the 
internal and external radii of the cylinders. The main 
conclusions of our work are addressed in section IV. 



II. THEORY 

Take a periodic array of infinitely long multishell nan- 
otubes of inner and outer radii r and R, respectively, ar- 
ranged in a square array with lattice constant a = 2 x R, 
as shown in Fig. 1. These tubules are assumed to be 
embedded in an insulating medium, with a real and pos- 
itive dielectric constant e$. In the energy range of inter- 
est in the interpretation of absorption cross sections and 
energy-loss spectra the diameter of typical multishell car- 
bon nanotubes (2 R ~ 10 nm) is small in comparison to 
the wavelength of light, and we also assume that this 
diameter is large enough that a macroscopic dielectric 
function is ascribable to the tubules. For simplicity, the 
magnetic permeabilities will be assumed to be equal to 
unity in all media. 

Planar graphite is a highly anisotropic material, and 
the dielectric function is a tensor. This tensor may be di- 
agonalized, by choosing cartesian coordinates with two of 
the axes lying in the basal plane and the third axis being 
the so-called c-axis. One defines the dielectric function 
£±(u>) perpendicular to the c-axis and the dielectric func- 
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tion £||(w) for the electric field parallel to the c-axis. For 
carbon nanotubes, we assume full transferability of the 
dielectric tensor of planar graphite to the curved geome- 
try of carbon tubules, as suggested by Lucas et aZEH for 
the case of multishell fullerenes. Hence, we simply assign 
a local graphite- like dielectric tensor to every point inside 
the nanotube and outside the inner core, and write 



e(uj) = e ± (uj){99- 



(w)rr, 



(1) 



where 99, zz, and rr are the unitary basis vectors of 
cylindrical coordinates. 

In the long- wavelength limit, a composite material may 
be treated as if it were homogeneous, with the use of an 
effective dielectric function £ e //. The optical absorption 
cross section of the composite is then directly given by 
lme e f f(u>). Also, for small values of the dimensionless 
parameter qR (qR < 1), q being the momentum transfer, 
the energy-loss spectra of a broad beam -of. swift elec- 
trons penetrating the composite is foundlifrEa to be well 
described by the q — > limit of the imaginary part of the 
effective dielectric function, i.e., the so-called energy-loss 
function, Im[—e~^(u))]. 

We consider an electromagnetic wave normally inci- 
dent on the structure, so that k y = k z = 0. For this 
propagation direction there are two different values of 
e e ff(ui) corresponding to s and p polarizations. In the 
case of s polarization the electric field is parallel to the 
cylinders at every point, and is not modified by the pres- 
ence of the interfaces. Hence, the s effective dielectric 
function of the composite is simply the weighted average 
of the dielectric functions of the constituents.E3 

For electromagnetic waves polarized normal to the 
cylinders (p polarization), the electric field may be 
strongly modified by the presence of the interfaces. An 
elementary analysis shows that 

( £e// -£o)E = /(£-£ )-E m , (2) 

where E is the average electric field over the composite, 

E = /E in + (1 - f)E out , (3) 

E in and E out representing the average electric field inside 
and outside the tubules, respectively, both lying in the 
plane of periodicity. 

In the case of a single two-dimensional circular inclu- 
sion (plain or hollow cylinder) embedded in an otherwise 
homogeneous medium (/ — > 0), one easily finds 



£o«E, 



(e-eo) -Ej 
and Eq. ®) then yields 

£ e ff = £o (1 + fa) 



(4) 



(5) 



where a represents the in-plane dipole polarizability per 
unit volume. 

As long as the composite is made of a sparse (/ << 
1) distribution of cylinders the presence of multipolar 



modes can be neg lectedji and the interaction between 
the cylinders can be introduced by simply replacing the 
average electric field E in Eq. (||) by E out . Then, with 
the aid of Eqs. (||) and (||) , one finds 



= E0 1 + / 



1-fLa 



and 



£ eff ~ £ 



1-/ 



1 + fLa 



(6) 



(7) 



with the .geometrical factor L = 1/2. For cylindrical 
geometry,E3 



(e,|A - £„)(£,, A + e ) - (e,, A - e )(ei|A + e )p 2A 



(£|| A + £ )(£|| A + £ ) - (£|| A - £ )(£|| A - E )p 

where 



2A ■ 



(8) 



(9) 



Eq. (g) with p = and s\\(w) = £±(<j-0 reduces to 
the well-known polarizability per unit volume of a plain 
isotropic cylinder. 

Eq. (^) is a generalization, appropriate for anisotropic 
tubules, of the well-known MG effective dielectric func- 
tion deri-ued by Maxwell-Garnett for a system of spherical 
particlcs.tj 

In the non-sparse or packed regime, where the presence 
of multipolar modes cannot be neglected, the inclusion of 
the full electromagnetic interaction between the tubules 
is unavoidable. In order to compute, with full inclusion 
of this interaction, the effective dielectric function of our 
periodic svstem, we have followed the method developed 
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and |24| for the calculation of dispersion re- 
lationships k(u>) of Bloch waves in structured materials 
with tensor-like dielectric functions. 

In the long-wavelength limit the composite material 
supports, for each polarization, only two degenerate elec- 
tromagnetic Bloch waves with vectors k and — k for which 
k(u>) roughly follows the dispersion relation of free light. 
Hence, in this limit the composite may be treated as if it 
were homogeneous, with an effective transverse dielectric 
function 



k 2 {oj)c 2 



(10) 



where c represents the speed of light. 

III. RESULTS AND DISCUSSION 

We consider a periodic array of hollow multishell car- 
bon nanotubes (see Fig. 1) with £ = 1, and take the 
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principal dielectric functions £j_(w) and £|i(w) of graphite 
from Ref . |25|. The real and imaginary parts of these com- 
ponents, as well as their corresponding energy-loss func- 
tions, are represented in Fig. 2 for energies up to 8eV. 
In this energy range, optical transitions mainly involve 
the 7r bands arising from the atomic 2p z orbitals. The 
peak in Im£j_(o;) at ~ 4.6 eV is associated with the max- 
imum in the joint density of states of the tt valence and 
conduction bands. The so-called tt plasmon at ~ 7eV, 
where Kee±(uj) and bae^(w) are small, is due to a tt—tt* 
interband transitionrS 

For energies over ~ 5 eV, the principal dielectric func- 
tion £j_(w) is close to that of a free-electron gas. Thus, 
for these energies tt electrons within each shell act as if 
they were free. However, for free-electron-like materials 
Kce±(u>) approaches large negative values with decreas- 
ing energy, whereas for graphite the plasmon region with 
Re e±(u>) negative is bounded from below and from above 
at about 5 and 7eV. 

Calculations of the s component of the effective di- 
electric function of an array of coaxial nanotubes were 
reported in Ref. [l5| for various values of the filling frac- 
tion / of the tubes, showing that they roughly reproduce 
the experimentally determined Ime e //(w) for / ~ 0.5. 

Here we focus on the case of electromagnetic waves 
polarized normal to the cylinders (p polarization) , which 
we first investigate within the MG approach. 



A. Maxwell-Garnett approach 

First of all, we ignore the anisotropy and simply as- 
sume that e(uj) = £\\(u>) — £±(lj). In this case, the in- 
plane dipole polarizability of Eq. (||) can be expressed in 
the form of a spectral representation: 



(11) 



it — m_ u — m+ 1 
where u is the spectral variable 

U= (l-£/£ r\ (12) 

m± are depolarization factors, 

1 



m± = - (1 ± p) 



(13) 



and B± are the strengths of the corresponding normal 
modes, 
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(14) 



The depolarization factors of Eq. (|13j) lie on the segment 
[0,1] for all values of p. Hence, normal modes only occur 
in the so-called plasmon region where the dielectric func- 
tion s{uj) is negative or equal to zero. For graphite, this 
plasmon region is bounded from below and from above at 



about 5 and 7eV where the spectral variable u is never 
smaller than u ~ 0.2. As the ratio between the inner and 
outer radii of the tubes goes to unity (p — > 1), one easily 
finds 



(e - e -i) 



(15) 



In the case of an isolated plain (p = 0) or hollow 
(p 7^ 0) cylinder, both Ime e //(o;) and Im[— e~^(uj)] are 
proportional to the imaginary part of the in-plane dipole 
polarizability a, and satisfy the relation 



Irn 



— £ 



eff 



= e n Ime 



eff- 



(16) 



Hence, the absorption and energy-loss spectra of iso- 
lated plain (p = 0) cylinders in vacuum (eo = 1) ex- 
hibit a strong maximum at ~ 6.5 eV where u = 1/2 
and e(ui) = —1, as shown in Fig. 3. For hollow tubes 
(p 0), there are two distinct dipolar modes with cither 
tangential or radial symmetry, at u = rn_ < 1/2 and 
u = m + > 1/2, respectively, similar to those present m 
the case of a thin planar filmEZl and a spherical shcll.Ej 
The corresponding peak heights in Im a(u) are easily 
found to be H B± j y/m± , H representing the peak height 
in the bulk energy-loss function, Im[— e _1 (w)]. Thus, the 
tangential mode at u = (1 — p)/2 appears to be more pro- 
nounced than the radial mode at u = (1 + p)/2 (see Fig. 
3). For p = 0.2 these modes are not resolved, due to the 
smoothing effect of the large damping associated with 
non-negligible values of lme±(uj), and the effect of the 
empty core is simply to redshift and soften this combined 
plasmon mode. For values of p in the range 0.2 < p < 0.8 
the tangential and radial plasmons are clearly identified. 
However, for larger values of p the tangential resonance 
condition (u — > 0) is never satisfied. Furthermore, in the 
limit as p —> 1 the effective dielectric function is given by 
Eq. (|l5|), thus showing the radial plasmon resonance at 
~ 7eV, where u = 1 and e(u>) — 0, and the character- 
istic peak at ~ 4.6 eV associated with the maximum in 
Ime(cj). 

Still ignoring the anisotropy, the MG effective dielec- 
tric function and inverse dielectric function of Eqs. (|^) 
and (ffl) can also be expressed in the form of a spectral 
representation: 



?eff - £o 



and 



e eff - £ o 



1-/ 



1 + / 



B, 



B- 



m,- 



B. 



B 4 



u — n_ u — n + 



where 



m ± = Ul- l -f± l -JpTXp 1 



n± 



1 



4p< 



(17) 



(18) 



(19) 



(20) 
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and 



B+ = - 



1 y // 2 + V ± f 
VP + V ' 



(21) 



the mode strengths (_B_ < 1/2 and £?+ > 1/2) adding 
up to unity, i.e., B+ + £>_ = 1. Also, m_, n_ < 1/2 and 
m+, n + > 1/2. 

In Fig. 4, we show the surface-mode positions m± and 
n± of Eqs. ( |l9| ) and (p0|), as a function of p, for various 
values of the ratio x between the lattice constant and the 
outer diameter of the cylinders. In the dilute limit (x — > 
oo), both m± and n± coincide with the depolarization 
factors of Eq. (|l^) entering the spectral representation 
of the polarizability. As x decreases, for each value of p 
these surface modes split into four distinct modes, m± 
and n±, which satisfy the simple relation: 



n± = 1 



(22) 



The strengths B± , as obtained from Eq. (|21j) , are plot- 
ted in Fig. 5, as a function of p, also for various values 
of the ratio x between the lattice constant and the outer 
diameter of the cylinders. As in the case of the depolar- 
ization factors, in the dilute limit (x — > oo) the strengths 
B± coincide with the strengths B± = 1/2 entering the 
spectral representation of the polarizability. As x de- 
creases, these strengths are different, especially for small 
values of p, showing that the strongest modes occur at 
u = rri- and u — n + with the depolarization factors m_ 
and n+ lying on the segments [0,(l-f)/2] (low-energy res- 
onance) and [(l+f)/2,l] (high-energy resonance) , respec- 
tively (see Fig. 4). The modes at u = m,- and u = n +1 
which have no strength for p = 0, lie on the segments 
[1/2,1] and [0,1/2], respectively. For x — ► 1, one easily 
finds B+ = 1/(1 + p 2 ) and B_ = p 2 /{l + p 2 ), thereby 
showing that modes occurring at u = rti- and u = n + 
dominate. Hence, when cylinders are touching and for 
most values of p there is a single dominating resonance 
condition at u — ► (low energy) and u — > 1 (high en- 
ergy) in Ime e //(u) and Im[— £~ff{w)], respectively, the 
strength of the remaining resonance being negligible. 

From Eqs. ( |l7|) and (^8|), we have calculated 
Im£ e //(w)// and ha[—e~^(u))]/f for various values of 
x and the ratio p between the inner and outer radii of 
the cylinders. While for x — 3 both Im£ e //(w)// and 
Im[— ejjy(a;)]// are found to nearly coincide with the 
imaginary part of the polarizability of Eq. ( [Tl] ) (see 
Fig. 3), we have found that the trend with decreas- 
ing the distance between the cylinders is for the low- 
energy and high-energy resonances to slightly move from 
the single-cylinder dipole resonances at u — (1 ± p)/2 
to lower and higher energies, respectively; this is obvi- 
ous in Fig. 6, where lme e f f(uj)/ f and Im[— e~ ff (u))]/f, 



as obtained from Eqs. (17) and (18), have been repre- 
sented for x = 1.5 and various values of p. We also note 
that as x decreases the low-energy (u = rn_) and high- 
energy (u = n+) modes, whose energy position depends 



only weakly on p, dominate the optical absorption and 
energy-loss, respectively. 

The role of anisotropy is displayed in Figs 7 and 8. 
In these figures, we have plotted Ima(w) (Fig. 7), 
lme eff (u)/f (Fig. 8), and Im[-e^(w)]// (Fig. 8), as 
obtained from Eqs. (||), (||), and (|7j) with full inclusion 
of the anisotropic dielectric function of graphite and for 
various values of p, as in Figs. 3 and 6. With the pres- 
ence of anisotropy, the two-mode structure exhibited by 
Eqs. (|ll]), (|l7|), and ( |l8| ) is replaced by a more compli- 
cated spectral representation. One sees from Fig. 7 that 
the tangential plasmon peak-position of isolated cylin- 
ders remains fairly insensitive to the anisotropy of car- 
bon nanotubes. However, its shape drastically changes, 
as a consequence of the presence of a nearly constant and 
positive dielectric function en, and the radial resonance 
condition is not visible. For x — 3 the MG dielectric func- 
tion is found to nearly coincide with the isolated-cylinder 
result, as in the isotropic case. As the distance between 
the cylinders decreases, the impact of the anisotropy is 
still to soften the resonances. Furthermore, one sees from 
Fig. 8 that now at small values of p the low-energy reso- 
nance is only visible in the optical absorption, while the 
energy-loss only exhibits the high-energy peak. 



1. Packed regime 

In this section we present results of our calculation 
of the effective dielectric function of a periodic array 
of carbon nanotubes, as obtained from Eq. ( |To| ) with 
full inclusion of both the anisotropy and the electromag- 
netic interaction between the tubules. All calculations 
presented here have been found to be insensitive to the 
precise value of the number of mesh points in the unit 
cell. For metallic structures, sampling meshes as large 
as 180 x 180 have beep_found to be required to provide 
well-converged results;E3 however, for carbon nanotubes 
sampling meshes of 60 x 60 have been found to provide 
well-converged results, which is due to the smoothing ef- 
fect of the large damping originated with the presence of 
interband transitions in graphite. 

In Figs. 9 and 10 we show our full calculations of the p 
component of the effective dielectric function of an array 
of plain (p — 0) carbon nanotubes, as obtained for vari- 
ous values of the ratio x between the lattice constant and 
the outer diameter of the cylinders: x = 2.0, x = 1.5 and 
x = 1.3 in Fig. 9, and x = 1.1 and x = 1.03 in Fig. 10. 
Also plotted in these figures by dotted lines are calcula- 
tions of the MG dielectric function of Eqs. @ and (0). 
These calculations show that the actual effective dielec- 
tric function is well described by the MG approximation 
in the low-filling-ratio regime, our full calculations be- 
ginning to deviate from those obtained within the MG 
approximation at x ~ 1.3. For smaller concentrations 
of graphite, multipolar modes cannot be neglected and 
the dipole-resonance positions necessarily deviate, as dis- 
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cussed in Ref. [l9]. 

One sees from Figs. 9 and 10 that, as within the MG 
approach, the trend with increasing the concentration 
of tubules is for the actual dipolar peak in Ini£ e ^(w) 
and Im[—e~ff(u))] to be shifted from the isolated-cylinder 
dipole mode at ~ 6.5 eV to lower and higher energies, 
respectively. When the nanoparticles are brought into 
close contact, electromagnetic coupling between them 
converts the dipolar surface mode into a very localized 
one, trapped in the region between the nanostructures, 
and the MG approximation fails to describe the details 
of the effective dielectric function. We note that multi- 
polar resonances, which are present in the close-contact 
regime, are not visible in the spectra due to the smooth- 
ing effect of the large damping characteristic of graphite. 
At higher concentrations of tubules, when graphite forms 
a connected medium (x < 1), dipolar modes cannot be 
excited. Hence, optical absorption exhibits a single peak 
originated in the maximum of Ime(o;) at ~ 4.6 eV, and 
the energy-loss function shows a single peak at the bulk 
plasmon resonance at ~ 7eV. Besides these peaks, there 
is a background of unresolved multipole contributions to 
Im £ e ff(u>) and hn.[—e~^(u))] in the plasmon region with 
energies between ~ 5 and ~ 7eV. 

Calculations of the p component (p polarization) of 
the effective dielectric function of a periodic array of hol- 
low (p 0) carbon nanotubes were reported in Ref. [l5| 
for various ratios p between inner and outer radii of the 
tubules. As in the case of plain cylinders (p = 0), these 
calculations nearly coincide for x = 2.0 with the results 
obtained within the MG approximation. However, in the 
close-packed regime (x = 1.03) the strong electromag- 
netic coupling between the tubes results in non-negligible 
contributions from multipolar resonances. This multipo- 
lar coupling provokes a redshift and a blueshift of the MG 
dipolar resonances that are visible in the optical spec- 
tra and the energy loss, i.e., the low-energy (u < 1/2) 
dipolar mode with tangential symmetry and the high- 
energy (u > 1/2) dipolar mode with radial symmetry. 
A comparison between these calculations and the exper- 
imentally determined macroscopic-p dielectric function 
of close-packed carbon nanotubea!3 was also presented 
in Ref. [L5|, showing an excellent agreement for p = 0.6, 
which yields in the close-packed regime (x = 1.03) a fill- 
ing fraction / ~ 0.5. 

Finally, we look at the E-field pattern associated with 
the electromagnetic resonances that are present in the op- 
tical absorption and the energy loss. In Fig. 11, we show 
detailed pictures of the intensity of the electric field and 
the corresponding charge density generated by normally 
incident p-polarized light impinging on our periodic ar- 
ray of hollow (p = 0.6) carbon nanotubes with x = 2. 
The frequency of the incident radiation has been chosen 
to be u)\ — 5.2 eV and 0J2 = 5.9 eV, corresponding to the 
location of resonances in the optical absorption and the 
energy loss, respectively. In this figure one clearly sees 
that the incident radiation is exciting prototypical dipo- 



lar tangential and radial dipolar plasmons at the surfaces 
of a nearly isolated hollow cylinder. At u>\ = 5.2 eV, the 
intensity of the electric field is maximum outside the nan- 
otube and the induced charge is clearly located on the 
outer surface of the tubule, showing a tangential field 
pattern. At 102 — 5.9 eV, the E-field intensity is maxi- 
mum in the hollow core of the nanotube and the induced 
charge is located on both surfaces of each tubule, though 
the largest part of it is located on the inner surface of 
the tube. These represent the main characteristics of a 
radial dipolar mode. 

In the case of a close-packed structure, the incident 
radiation excites a very localized tangential surface plas- 
mon. This is illustrated in Fig. 12, where we show 
pictures of the E-field intensity and the correspond- 
ing charge density generated by normally incident p- 
polarized light impinging on our periodic array of hol- 
low (p — 0.6) carbon nanotubes with x = 1.03. At 
lui = 4.75 eV, where the optical absorption exhibits a 
maximum, the E-field intensity is strongly enhanced in 
the region between the tubes and the induced charge is 
clearly localized, showing a strongly localized tangential 
field pattern. At u>2 = 6.1 eV, where the energy loss is 
maximum, one clearly sees that the incident radiation is 
exciting a radial surface plasmon which is very similar 
to that observed in the case of nearly isolated nanotubes 
(see Fig. 11). 

IV. SUMMARY AND CONCLUSIONS 

We have reported extensive calculations of the effec- 
tive electronic response of aligned multishell carbon nan- 
otubes, of interest in the interpretation of absorption 
spectra and electron energy-loss experiments. A local 
graphite-like dielectric tensor has been assigned to every 
point of the multishell tubules, and the effective dielectric 
function of the composite has been computed by solving 
Maxwell's equations in media with tensor-like dielectric 
functions. A MG-like approach appropriate to the case 
of infinitely long anisotropic tubules has also been devel- 
oped, showing that MG results are accurate as long as the 
distance between the axis of neighboring plain tubules is 
not smaller than ~ 1.3 times the outer diameter of the 
cylinders. 

The effective response of carbon nanotubes has been 
found to be sensitive to both the inner cavity of hollow 
tubules and the anisotropy, which we have first inves- 
tigated with use of the MG approximation. Within this 
approach, we have analyzed the mode strengths and posi- 
tions of tangential and radial surface plasmons, and have 
investigated the effect of the anisotropy and the electro- 
magnetic interactions between the tubes. 

Finally, we have presented an analysis of the electric- 
field pattern associated with the electromagnetic reso- 
nances that are present in the optical absorption and the 
energy loss. We have shown that incident p-polarized 
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light impinging on a periodic array of hollow carbon nan- 
otubes excites tangential and radial plasmons at the sur- 
faces of each tubule. In the case of isolated hollow nan- 
otubes, these surface plasmons are found to be prototyp- 
ical dipolar modes. In the close-packed regime, where 
nanotubes are nearly touching, our calculations indicate 
that the incident radiation excites a very localized tan- 
gential surface plasmon, while the radial plasmon is found 
to be very similar to that observed in the case of isolated 
nanotubes. 
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FIG. 1. Multishell nanotubes of inner and outer radii r and 
R, respectively, arranged in a square array with lattice con- 
stant a. The cylinders are infinitely long in the y direction. 
The electromagnetic interaction of this structure with a nor- 
mally incident plane wave of momentum k [k y = k z — 0] and 
energy u is investigated. 

FIG. 2. Energy dependence of the real and imaginary parts 
of the principal dielectric functions, Ree(o;) and Ime(w), and 
the energy-loss function, Im[— e _1 (w)], for graphite, taken 
from Rcf. ^ Solid and dashed lines represent components 
perpendicular and parallel to the c-axis, respectively. 

FIG. 3. Imaginary part of the in-plane dipole polarizability 
per unit volume of isotropic plain (p = 0) and hollow (p 0) 
cylinders, as obtained either from Eq. (^) with e\\ (uj) — e±(uj) 
or from Eq. ( J5~T| ) , as a function of frequency lu and for various 
values of the ratio between the inner and outer radii of the 
cylinders: p = 0, 0.2, 0.4, 0.6, and 0.8. 

FIG. 4. Surface-mode positions: (a) m- ,n+ and (b) 
m+, n_ , as obtained from Eqs. ([H)|) and (|5^) versus p, for var- 
ious values of the ratio between the lattice constant and the 
outer diameter of the cylinders: x = 3 (solid lines), x = 1.5 
(dotted lines), and x = 1 (dashed lines). 

FIG. 5. Strengths B±, as obtained from Eq. (^l|) versus 
p, for various values of the ratio between the lattice constant 
and the outer diameter of the cylinders: x = 3 (solid lines), 
x — 1.5 (dotted lines), and x = 1 (dashed lines). 
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FIG. 6. Imaginary part of the MG effective dielectric func- 
tion and MG effective energy-loss function of the periodic 
system described in Fig. 1, as obtained from either Ens. ([]) 
and (fit) with e\\(uj) = £±(u) or from Eqs. (|l|) and Q, and 
for various values of the ratio between the inner and outer 
radii of the cylinders: p — (solid line), p = 0.2 (dotted line), 
p — 0.4 (short-dashed line), p — 0.6 (long-dashed line), and 
p — 0.8 (solid line). The ratio between the lattice constant 
and the outer diameter of the cylinders is taken to be x — 1.5. 

FIG. 7. Same as Fig. 3, as obtained from Eq. (^) with use 
of the actual principal dielectric functions e±(u>) and eii (u) of 
graphite. 

FIG. 8. Same as Fig. 6, as obtained from Eqs. Q and (j?|) 
with use of the actual principal dielectric functions e±(u>) and 
£|| (w) of graphite. 

FIG. 9. The real and imaginary parts of 
the long- wavelength effective dielectric function, Ree e //(w) 
and lxas e ff(uj), and the energy-loss function, lm[—e~}Ju)], 
of a periodic array of plain (p = 0) carbon nanotubes, for 
p polarized electromagnetic excitations. Solid, long-dashed, 
and short-dashed lines represent our full calculated results, 
as obtained from Eq. ([to]) for ratios between the lattice con- 
stant and the outer diameter of the cylinders x = 2, x = 1.5, 
and x — 1.3, respectively. The dotted lines represent MG 
results, as obtained from Eqs. ([]) and ([?j). 

FIG. 10. Same as Fig. 9, for ratios between the lattice 
constant and the outer diameter of the cylinders x — 1.1 
(solid line) and x = 1.03 (dashed line). 

FIG. 11. The intensity of the electric field and the cor- 
responding charge density generated by a normally incident 
p polarized electromagnetic plane wave impinging on a pe- 
riodic system of hollow (p = 0.6) carbon nanotubes with 
x — 2.0. Both the electric field and the induced charge have 
been evaluated at two different frequencies: Wi = 5.2 eV (the 
frequency at which the optical absorption is maximum) and 
UJ2 = 5.9 eV (the frequency at which the energy loss is maxi- 
mum). White/dark areas mean small/large values of the elec- 
tric-field intensity and the induced charge. 

FIG. 12. Same as Fig. 11, for x = 1.03, wi = 4.75 eV, and 
uj2 — 6.1 eV. 
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